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Abstract 

We systematically investigate the mode dispersion and spectral weight of 
the elementary excitation spectra in one-dimensional continuum and lattice 
electron systems by using the RPA, the Luttinger liquid model, and the Hub- 
bard model. Both charge and spin excitations are studied in details and 
compared among the theoretical models. For the lattice Hubbard model we 
use both Bethe-ansatz equations and Lanczos-Gagliano method to calculate 
dispersion and spectral weight separately. We discuss the theoretically cal- 
culated elementary excitation spectra in terms of the experimental inelastic 
light (Raman) scattering spectroscopy of ID semiconductor quantum wire 
systems. Our results show that in the polarized (i.e. non-spin-flip) Raman 
scattering spectroscopy, only the ID charge density excitations should show 
up with observable spectral weight with the single particle excitations (in 
RPA) or singlet spin excitations (in the Luttinger model and the Hubbard 
model) having negligible spectral weight. The depolarized (spin-flip) Raman 
scattering spectra manifest the spin density or the triplet spin excitations. We 
also provide a qualitative comparison between the continuum and the lattice 



ID elementary excitation spectra. 
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I. INTRODUCTION 



There has been considerable recent interest in the elementary excitation spectra of one- 
dimensional (ID) systems, both in continuum systems such as ID semiconductor quantum 
wires (QWR) or carbon nanotubes and in lattice models such as ID organic conductors. 
In some sense the ID elementary excitation spectrum is remarkably simple |T]-§| because 
single-particle (electron-hole like) excitations are completely suppressed in one dimension, 
and only ID collective (charge density and spin density) excitations exist at long wavelengths. 
In practice, however, resonant inelastic light (Raman) scattering (RRS) based spectroscopic 
studies of ID semiconductor (GaAs) quantum wires find several electronic excitation 
modes in polarized and depolarized Raman spectra whose mode dispersions and spectral 
weights have been measured carefully. In the polarized geometry incident and scattered pho- 
tons have the same polarization vectors indicating the absence of spin flips in the electronic 
excitations, while in the depolarized geometry the incident photons change the polarization 
after scattering due to the spin flips of electrons. From the many-body linear response 
theory, the polarized spectrum is proportional to the imaginary part of the charge density 
correlation function and the depolarized spectrum is proportional to the imaginary part of 
the spin density correlation function. One particularly noteworthy feature of these Raman 
scattering studies has been the finding that the qualitative features of ID GaAs QWR ele- 
mentary excitation spectra (in both polarized and depolarized geometries) are remarkably 
similar to those in 2- and 3-D GaAs systems. 

The goal of this paper is to investigate theoretically the charge and spin elementary exci- 
tation spectra as well as their spectral weights in ID electron systems, both for a continuum 
jellium electron gas and for an atomic lattice model. Our calculations (in particular, our 
spectral weight calculations) should apply directly to the experimental Raman scattering 
data if the resonance effects are dynamically unimportant in the interpretation of the Ra- 
man experiments. We refer to our calculations as the nonresonant Raman scattering (NRS) 
theory where only the conduction band electrons are taken into account as opposed to the res- 
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onant Raman scattering where both the conduction band and the valence band participate. 
In fact, the theory developed in this paper has been the standard theory for discussing the 
resonant Raman scattering spectroscopy until very recently when several publications f^HH 
dealing with the full subtleties of the resonance effect have appeared in the literature. We 
emphasize that, quite apart from the resonant Raman scattering spectroscopy, theoretical 
results presented in this paper stand on its own as a comprehensive theory for the elementary 
excitation spectra of ID electron systems. 

In this paper, we will study the standard (nonresonant) Raman scattering spectroscopy 
in three theoretical models: the RPA-Fermi liquid model, the Luttinger liquid (LL) model, 
and the ID Hubbard model. As emphasized above, by "nonresonant" we mean that the 
theory neglects all effects of the valence band in resonant Raman scattering (which is a 
2-step process, with the incident photon being absorbed by a valence band electron which 
thereby gets excited into an excited conduction band state with an electron from inside 
the conduction band Fermi surface subsequently combining with the valence band hole 
with the emission of the scattered photon). If the valence band can be ignored, then only 
conduction band density fluctuations are responsible in the linear response theory of the 
scattering process. The calculation is then simplified to the evaluation of the density- density 
correlation function (for polarized spectrum) and the spin-density correlation function (for 
depolarized spectrum), whose imaginary parts are proportional to the spectra measured in 
the experiments. This approach of identifying the measured elementary excitation spectra in 
the Raman scattering experiments as the charge (polarized spectra) or the spin (depolarized) 
density correlation function of the electron system in the conduction band has a long and 
fairly successful history |H| in the semiconductor structures. We take the same approach 



here, and construct our charge and spin density correlation functions (which give the spectral 
strengths of the elementary excitations through their imaginary parts or the corresponding 
dynamical structure factors) entirely from the conduction band carriers, ignoring all effects 
of the valence band in the resonance process. 

There has been one persistent feature in the experimental Raman spectra of semicon- 



ductor systems, including ID QWR structures, which does not seem to have an obvious 
explanation in terms of the non-resonant theory discussed in this paper. There is often a 
low energy spectral peak in the polarized spectra at an energy well below the expected collec- 
tive charge density excitation (CDE) peak (and in addition to the charge density excitation 
peak, which always shows up at the usual energy). This additional peak occurs around the 
single particle excitation energy, which typically contributes little to the dynamical struc- 
ture factor (i.e. the density correlation function) at the low wavevectors (q <C fop) of Raman 
scattering experiment, and therefore should have negligible (unobservable) spectral weight. 
There have been many suggestions for the resolution of this puzzle (namely, why the sin- 
gle particle excitation weight is enhanced in the density correlation spectrum), and we will 
quantitatively consider several of these suggestions in this paper. Our conclusion, based on 
the results presented in this paper, is that this puzzle in all likelihood arises necessarily from 
the resonant nature of Raman scattering experiments, as has recently been argued in the 
literature which is beyond the scope of this paper. Our critical quantitative consider- 

ation of the several suggested scenarios (within the non-resonant theory of using conduction 
band properties only) for explaining why the single particle excitation has large spectral 
weight shows that none of them is capable of resolving this problem quantitatively. While 
studying the ID elementary excitation spectra is the primary goal of this paper, considering 
the single particle excitation spectral weight issue in the non-resonant Raman scattering is 
one of our important secondary goals. 

We first present the results of the Fermi-liquid random phase approximation (RPA) 
calculation for this problem in Sec. II. The RPA calculation has been shown to give a good 
description for the dispersion relations of the elementary excitations in comparison with 



the experimental results for both intersubband and intrasubband ID excitations ||12|| . 
Being a standard Fermi liquid (FL) theory, however, the RPA calculation is unable to explain 
the relatively large spectral weight of the "single particle excitation" (SPE) in the polarized 
spectrum of the experiment as discussed above. We include the effects of the breakdown of 
momentum conservation and the nonparabolic energy dispersion in our calculation in Sec. 



II-D to check if they can explain the SPE feature within RPA, but neither gives qualitatively 
correct results for the polarized spectrum. In the Luttinger liquid model we present in Sec. 
Ill, we find zero weight at the SPE energy (as we should, since in the LL model there are no 
single particle like quasi-particle excitations), and all the spectral strength is at the charge 
boson mode, which is exactly the CDE mode in RPA. In Sec. IV, we use the ID (lattice) 
Hubbard model with repulsive on-site spin- dependent interaction to study this problem. 
The study of ID elementary excitations and the associated spectral weights in the Hubbard 
model is one of our main results in this paper. This Hubbard model study was originally 
motivated by the suggestion in ref. that a possible way to interpret this SPE puzzle 
(the existence of a single particle peak in Raman scattering) in ID should be different from 
those in higher dimensions, and the so-called SPE peak may be arising from the spin-singlet 
excitations (SSE) of interacting ID systems |r3|JT^] . Therefore, we choose to study in detail 
the ID Hubbard model, in which the spin- dependent interaction is expected to enhance the 
contribution of the spin-singlet excitation, which is proposed [ |T3| , |T4f as the extra SPE-like 
feature showing up in the experiment. Although the Hubbard model is a lattice model 
and consequently may not apply directly to the continuum QWR system, we argue that 
it is useful to understand the detailed excitation spectra in the ID Hubbard model in the 
context of this problem because one can quantitatively study the interacting ID elementary 
excitation spectra using the Hubbard model. 

Our conclusion based on our work (as presented in Sec. II-IV) is that all of these NRS 
theories fail in providing an explanation for the large "SPE" spectral weight in the long 
wavelength region observed in the experimental polarized RRS spectrum, although there 
could be substantial spectral weight in the single particle excitation at large wavevectors 
(short wavelengths). The spectra of depolarized nonresonant Raman scattering is also stud- 
ied in these different models, and good general agreement with experiment is found in the 
sense that the depolarized spectra describe the spin density excitation (SDE) or the triplet 
spin excitation. 
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II. FERMI LIQUID MODEL 



In the standard many-body theory, which neglects all interband resonance effects and 
considers only the conduction band, the NRS intensity for polarized spectrum is proportional 
to the imaginary part of the charge density correlation function Xp- 

oc -Imx P {q,uj), (1) 



dVtduo 7r 

where 



oo 



X P (q,u)=i dte^ Xp (q,t), (2) 
Jo 

and 

Xp (q,t) = ([p\q,t),p(q,0)}) , (3) 

where p(q, t) = n(q, t)—n is the density fluctuation operator and n(q, t) = J2 ps c \>+q,s{t) c p,s{t) 
is the density operator (c PtS (t) is the electron annihilation operator). For the depolarized 
spectrum, we just need to change Xp t° Xa, the spin density correlation function, which is 
defined by 

Xa(q,t) = ([a\q,t),a(q,0)]) o , (4) 

where cr(q,t) = J2 ps s c p+q,s(t) c p,s(t) is the spin density operator. 

Since Xpiq^) is equivalent [T^] to the reducible polarizability, U p (q, u>), we may rewrite 



(Fig. 1) it as a sum of the series of irreducible polarizability, Uo tP (q, u>), bubble diagrams: 



n p (g, to) = n 0iP (g, u) + n 0iP (g, u) v c (g)n , P (g, w) + 

TL 0t p(q,u) Uo tP (q,u) 



(5) 



1 - V c (q)U 0iP (q,u) e(q,uj) 
where e(q, u) is the dynamic dielectric function. For the spin density correlation function, 

Xa{q, t), however, the equivalent reducible response function, H a (q, to) is exactly the same as 

the irreducible one, Ho t(T (q, to), due to the spin independence of Coulomb interaction (i.e. V c 

is explicitly spin independent). We will then use RPA [|nj to calculate H p (q, to) and Ilo, CT (g, to) 

and evaluate the charge and spin density correlation function (reducible polarizability) in 

the following sections. 



A. RPA and Hubbard approximation 



In the random phase approximation, which is very successful in the high density inter- 
acting two- or three-dimensional electron systems (the corresponding Feynman diagrams 
are shown in Fig. 1), the irreducible polarizability is approximated by the noninteracting 
electron-hole bubble (and therefore U 0jP / a (q,uj) = Itfj^fau)). We can calculate its ID 
analytical expression at T = (we choose h — 1 throughout) to be 



(27T) 2 

= — dp 

7T J 



m 
= — In 



irq 
m 
7iq 



dvdp Gq(p, v)Go(p + q,v + to) 

np(p) -n (p-q) 
ui + i5 — p 2 /2m + (p — q) 2 /2m 
(u + id) 2 - (q 2 /2m - qv F ) 2 ' 
(uj + i5) 2 - (q 2 /2m + qv F ) 2 _ 



In 



UJ 



(q 2 /2m — gyp)' 



(q 2 /2m + qvF)' 



+i sign(cjg) 



9 I — qvp — - — 



2m 



9 \ \ui\ — qvF + 



2m 



(6) 



where G (p, v) is the bare electron Green's function and n (p) = 9{kp — \p\) is the noninter- 
acting momentum distribution function. The first term of the last equality is the real part 
while the second term is the imaginary part of n^^.(g, uj), from which we could determine 
the single particle excitation region as shown in Fig. 2(a). We should note that unlike the 
2D or the 3D situation, the SPE continuum region for ID system is very restricted in the 
long wavelength limit (q k F ). In higher dimensions, the SPE continuum is gapless in 
the low energy region at any finite wavevector smaller than 2k F , but it is gapped in the 
ID system due to the phase space restriction. The strength of the SPE is a constant, inde- 
pendent of the frequency, uj, as we can see from Eq. (6), so that its contribution must be 
very limited in the interacting situation, which is dominated by the CDE or the collective 
plasmon excitations. From Eq. (6) we also see how the SPE region disappears completely 
as the band curvature goes to zero, i.e. when the band dispersion is taken to be linear as in 
the LL model. 



To improve the RPA calculation, which does not consider any vertex correction in 



H^p^ cr (q,uj), and therefore is exact only in the long wavelength limit, Hubbard intro- 
duced a simple local field correction to RPA by approximately summing the ladder exchange 
diagrams (and thus including vertex correction in a very crude manner). Instead of going 
through the details of his derivation, we here only show the final ID result, which is sim- 
ple |§ and can be applied to our calculation. The RPA dielectric function of Eq. (5) is 
replaced in the Hubbard approximation (HA) by 

^)-l- 1 + Vc{g)amzAM , (7) 

where the local field(vertex) correction, G(q), in the ID Hubbard approximation is given 
by I 



1 VJJq 2 + hi) 

In the long wavelength limit, q — > 0, G(q) — > for Coulomb interaction and the RPA result 
is restored in the long wavelength limit as it must. In Fig. 2 we plot the dispersion and 
spectrum of the ID charge density collective excitation (usually called the ID plasmon mode) 
within both the RPA and the Hubbard approximation — the plasmon mode is defined by 
the zero of the dielectric function and the intensity or the spectral weight is given by the 
imaginary part of the dielectric function (i.e. the dynamical structure factor JT5|]). 



B. Charge density correlation function — polarized spectrum 

Using the RPA result, Eqs. (5) and (6), we can obtain the density correlation function 
as well as the polarized NRS spectrum: 



lmIl p (q,Lj) 



ImH*f(^) 



1 - V c (q)^K F p A (Q^)\ + [V c (q)hnIl^ A (q,u 
Note that the imaginary part of is zero outside the SPE region. Therefore we expect 
a delta function like excitation at 1 — V c (q)KeIl^ ) A (q, uj) = for the plasmon mode or the 
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CDE dispersion. The dispersion is plotted in Fig. 2(a). At T = without any impurity 
scattering the dispersion of the plasmon mode in RPA is given by 



where A(q) = exp[qir /mV c (q)}. In the long wavelength limit (q — > 0), the dispersion becomes 



due to the fact that V c (q) — > ln(l/|g|) diverges in the long wavelength limit. Note that Eq. 
(11) is exactly the same as the dispersion obtained in the pure LL theory, which is discussed 
in Sec. Ill and is shown in Eq. (15) below. In Fig. 2(a), we find that the plasmon energy 
is actually larger than the SPE continuum energy for all momentum, so that there is no 
Landau damping in the ID system within the RPA calculation. The ID plasma dispersion 
has no gap in the long wavelength limit, but an infinite slope at q = due to the logarithmic 
divergence of the ID Coulomb interaction. 

In Fig. 2(b), we also show the typical calculated polarized RRS spectrum (proportional 
to Imfl in Eq. (9)) using a phenomenological broadening factor, 7 = O.ObEp, which may be 
arising from impurity scattering. We find that the spectral weight of CDE is much larger 
than that of SPE (about one thousand times !). In the same figure we show the HA results 
as well. We find that while the vertex correction indeed increases the SPE weight somewhat 
relative to the CDE weight, the HA is still completely unable (by a factor of 100) to explain 
the experimental finding ||] of the SPE mode being comparable in the intensity to the CDE 
mode (double-peak structure) in the polarized RRS spectrum [|J. 

Moreover, if the electron energy dispersion is linear (as it is close to the Fermi point), 
the SPE excitation spectral weight calculated by Eq. (9) disappears. This indicates that the 
band curvature around kp plays an important role in forming the SPE peak in the experi- 
ments. For example, in the linearized LL theory the SPE will have (Sec. Ill) exactly zero 
spectral weight. In Sec. II-D below we consider the band nonparabolicity effect explicitly. 



Vpi(q) 



A(q)(q 2 /2m + qvp) 2 — (q 2 /2m — qvp) 2 



(10) 




(ii) 
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C. Spin density correlation function — 



depolarized spectrum 



As discussed above, the depolarized spectrum of ID electron systems in the RPA is just 



maximum and minimum energies (for uj > 0), qvp + q 2 /2m and qvF — q 2 /2m respectively, 
the same as the SPE in the polarized spectrum. Thus, within RPA the depolarized mode 
dispersion is identical to the SPE energy, i.e. uj = qvp with a q 2 broadening. One should 
note, however, that when vertex corrections such as in the Hubbard approximation are taken 
into account, the spectra of the spin density excitation (SDE) will not be exactly the same 
as the SPE due to the vertex correction induced energy shift. In Fig. 3 we show the SDE 
spectrum obtained by calculating the spin density correlation function in RPA and also in 
the Hubbard approximation. The vertex correction shifts the SDE peak to lower energy (an 
excitonic shift) in the HA and thus separates it from the SPE mode. 

D. The breakdown of momentum conservation and nonparabolicity effects in RPA 

Beyond the standard RPA calculation, we include two nongeneric effects, the breakdown 
of momentum-conservation (arising from impurity scattering, for example) and the non- 
parabolicity of electron energy dispersion, because both of these corrections are likely to 
transfer some large wavevector SPE weight to smaller wavevector. For the breakdown of 
momentum-conservation, we use a phenomenological approach |T7[ by introducing a broad- 
ening function that couples the polarized spectrum x p (g, cu) at momentum q to that at 
momentum q'\ 



where T is a (phenomenological) dimensionless factor denoting the strength of the breakdown 
of momentum-conservation. For r — > we get back the original spectrum. In Fig. 4(a), 
we show the numerical calculation results of this effect by applying Eq. (12) onto the RPA 
result Eq. (9). At first sight, one finds that finite V does decrease the peak value of CDE 



the imaginary part of U^. (q, uj) in Eq. (6), which is nothing but a square function with 




(12) 
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and enhance the SPE weight. For T > 0.5, however, we find that the SPE peak merges 
into the very broad CDE peak, which is broadened also by the breakdown of momentum- 
conservation. In other words, the breakdown of momentum-conservation reduces the CDE 
peak strength and also broadens its width without changing either the total CDE spectral 
weight or the SPE weight qualitatively. Therefore, in our direct numerical calculation, we 
show that the breakdown of momentum-conservation (at least through phenomenological 
treatments, such as Eq. (12)) is not the candidate mechanism to provide an SPE spectral 
weight comparable to the CDE weight in RPA calculations. 

We now discuss the same issue by considering the band nonparabolicity effect of the 
electron energy dispersion. We recalculate the RPA spectral weight including band non- 
parabolicity via an additional q 4 term in the electron energy dispersion [18 



(13) 



This expression of E(q; A) keeps the Fermi energy constant (E(kp', A) = Ep) for all A and 
changes the electron effective mass m e (A) = m/(l — A) consistently. In Fig. 4(b), we show 
the calculated polarized RRS spectrum for different values of A < 0.1. We find that the 
enhancement of the SPE weight is very small, while the CDE peak almost keeps the same 
weight. Using larger A will cause greater blue-shifts in both SPE and CDE energies due to 
the increase of Fermi velocity via Eq. (13), which disagrees with the experimental results 
(note that the standard RPA results provide very good agreement with the experimental 
results in the excitation energies |4]]). Therefore, the nonparabolicity effect cannot enhance 
the SPE spectral weight to be comparable to the CDE weight. 



III. LUTTINGER LIQUID MODEL 

The Luttinger liquid model JTJfJ,|T^] is thought to provide a generic low energy description 
for ID electron systems, which are characterized by the LL fixed point in the renormalization 
group sense. The standard and exactly solvable LL model is the ID electron gas with a linear 
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dispersion (E(k) = rvp{k — rkp)) around Fermi points (±&_f) at each branch (r = ±1) and 
with short-ranged forward interaction |1],|2[. The Hamiltonian of this model is written as 

H = H v F (rp - k F ) : cl ps c rps : +# 4 + H 2 , (14) 

rps 

where : • • • : is the normal order operator, if 4 is for the electrons interacting within the right 
(r = +1) or the left (r = —1) band only and H2 is for the electrons interacting between 
these two bands. Using bosonization method and a linear transformation, the Hamiltonian 
of Eq. (14) can be exactly diagonalized by the two boson operators: charge boson p r (p) 
and spin boson cr r (p). This fact makes the collective excitation spectra (CDE and SDE) in 
the Luttinger model very simple: both the charge (CDE) and the spin (SDE) modes are 
delta function like poles and there is no SPE mode or equivalently, quasi-particle spectral 
weight, at all. In the following, we think it is instructive, for the sake of completeness and 
notational clarity, to present this result not only from the bosonization method we discussed 
above but also from the many-body Feynman diagrammatic technique via Ward identity, 
which is essentially equivalent to the RPA calculation but with linear band dispersion. 



A. Bosonization method 



Evaluating the expectation value, we get for the charge (p) and the spin (a) sectors: 



71 [U + lb) Z — {UJqf 



2L K a (q)v a (q)q 2 
Xv{q,u) = — 16 

where the excitation energy, uj^ u = \q\v p / a (q) = \q\v F / K p / a (q), and the charge sector Lut- 
tinger exponent K p (q) is 

Kp(q) J l + *YM\- l, \ (17) 

V TlVp J 

while K a (q) = 1 in the spin sector for the spin-independent Coulomb interaction. 

It is clear that the above results are completely spin-charge separated, which is another 
important feature of the LL model. One should note that there is no spectral weight in 
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Xp(Q, k>) at = g-yp for any SPE mode (or, for that matter any mode). This shows that 
the small SPE peak (compared to CDE) in the Fermi liquid-RPA theory is totally absent in 
the LL theory. Thus any possible explanation within the LL theory for the anomalous low 
energy peak in the polarized RRS spectra must arise from some mode (e.g. a multiboson 
mode or an SSE mode) other than the SPE mode which is completely absent in the LL 
theory. 

We note one other aspect (the spin-charge separation mentioned above) of the LL theory 
in this context which has created some minor confusion. The spin-charge separation of 
the LL theory has nothing whatsoever to do with the separate existence of SDE/CDE in 
the depolarized/polarized RRS spectra. The collective spin and charge density excitations 
are completely distinct excitations in the FL theory as well — they are the poles of the 
appropriate spin (for SDE) and charge (CDE) density correlation functions of the system 
which have totally different energies and selection rules (i.e. whether there is a spin-flip or 
not) in any reasonable theory. The reason spin-charge separation is rather complete in the 
LL theory is because the Luttinger liquid does not have any quasi-particles or single particle 
excitations — it has only collective spin and charge excitations which are poles of different 
correlation functions and are always separate. Indeed, higher dimensional systems, such as 
2D and 3D GaAs structures, exhibit qualitatively similar RRS spectra as in the ID system 
with the CDE peak (and a weak SPE-like low energy feature) showing up in the polarized 
spectra and the SDE peak showing up in the (spin-flip) depolarized RRS spectra. These 
higher dimensional systems are obviously Fermi liquids and have no LL-like intrinsic spin- 
charge separation while at the same time having distinct CDE and SDE collective modes. 

Unlike the formulae for the single-particle Green's function, in which the non-Fermi- 
liquid-like Luttinger liquid feature arises from the nonperturbative power-law behavior to- 
gether with the velocity renormalization, the charge and spin density correlation functions 
(which are two-particle Green's functions) have no such power-law behavior at all. The den- 
sity correlation functions and the the associated charge/spin density excitation collective 
mode spectra are essentially identical in the LL and the FL-RPA model H3"1,|TT| except for 
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the complete absence of any SPE spectral weight in the charge sector in the LL model. The 
Luttinger liquid effects appear only in the mode velocity renormalization, v p (q), and the 
overall mode amplitude factors K p (q) in the collective mode spectra. 



B. Diagrammatic method 

The diagrammatic method for the Luttinger liquid theory tells us more about the transi- 
tion from the Fermi liquid to the Luttinger liquid, because it is physically more transparent 
than the bosonization technique, which is more of a formal mathematical tool. Early sem- 



inal work by Dzyaloshinkii and Larkin |19| and recent important work of Schultz |20| have 
shown that this method is equivalent to the bosonization theory, even though its theoretical 
structure follows a Fermi liquid type conventional many-body theory. 

Eq. (5) is the starting point of the diagrammatic calculation. In order to evaluate the 
irreducible polarizability, we use the Ward identity connecting the Green's function with the 
vertex function in the following formula 

T rs (p, is, q, u) = — , (18) 

u> + 10 — rqvF 

where T rs (p, u, q, u) is the vertex function of two particle lines and one interaction line. The 
Ward identity follows directly from the particle and current conservation in each branch and 
spin (valid only for forward scattering) coupled with linear dispersion relation. It can be 
derived by summing the infinite series of vertex diagrams as shown in Fig. 5. Using this 
vertex expression, one can calculate the exact irreducible polarizability (consider the charge 
part only) 

U , p (q, u) = J 7^ J 7^ G rs(p, v)G rs (p -q,u-u) 



E 



xT rs (p, u,q,u) 
dpn rs {p-q) -n rs (p) 
2tt uj + iS — rqvF 
2q/n 



u + iS — rqvp 
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(19) 



Comparing Eq. (19) with Eq. (6) we find that they are identical if we only change the 
parabolic dispersion in Eq. (6) to the linear one as in the Luttinger model. Therefore we 
obtain the striking result that the irreducible polarizability of the linear band dispersion 
model is exactly the same as the RPA result. In other words, vertex corrections to the 
irreducible polarizability vanish. This result can also be verified by the topological argument 
given in ref. |]20[], which shows that all the electron- hole loops connecting with more than 
three interaction lines cancel out. Note that the above result is independent of temperature, 
due to the particle number conservation in the p-integral. 

Using this formula of Uo }P (q,u), the charge density correlation function can be easily 
obtained via Eq. (5). It is instructive, however, to see how it works in Eq. (9). In the usual 
RPA calculation in ID, the plasmon mode (given by 1 - U c (g)ReITg A = 0) is above the 
SPE region so that the square function (T = 0) of ImII^ A (g,w) still exists, unaffected by 
the CDE mode, while the CDE itself is a delta function like pole because lmIl^ A (q,u) = 
at the plasmon energy. When the dispersion is linearized as in the LL model, Imll^? (q, oS) 
itself becomes a delta function at u = qvp rather than a square function (see Eq. (19)), 
leading to the complete suppression and the disappearance of the SPE mode as we can 
see from Eq. (9) at uo = qvp. This makes this diagrammatic result consistent with the 
bosonization result, in which there is manifestly no single particle eigenstate at all in the 
final spectrum. Both approaches predict the complete absence of an SPE mode in the LL 
theory. 

Finally we can conclude that the LL theory gives precisely the same collective mode 
spectra as the corresponding RPA theory does with the only difference being the complete 
non-existence of any (even weak) SPE features in the LL spectra. 

IV. HUBBARD MODEL 

Motivated by the long-standing SPE-feature puzzle discussed in the Introduction, it has 



been suggested |I3],B| that we can interpret the "SPE" peak observed in the experimental 
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RRS spectra as the "singlet spin excitation" | 13f . In other words, the incident photon 



virtually flips the electron spin and then restores its polarization after the scattering, leaving 
the electron spin unchanged. Unlike the triplet spin excitation, SDE, which manifests itself in 
the depolarized RRS spectra, the virtual spin-flip process of SSE may, in principle, contribute 
to the final spinless scattering matrix element of the polarized Raman scattering spectrum, 
so one could expect that SSE should be very close to the SPE in energy under the spin- 
independent Coulomb interaction. That explains why one cannot simply separate these two 
(SSE and SPE) whether in the experimental measurement or in the theoretical calculation. 
In the LL theory, there is no SPE, but SSE is, in principle, allowed and may simulate the 
SPE of the FL-RPA theory. To investigate the role of SSE in more details we could use 
the ID Hubbard model on a lattice, which can be mapped to the exactly solvable Luttinger 
liquid model in the long wavelength limit. Thus the ID Hubbard model has no generic 
SPE properties, and could therefore be useful in the understanding of SSE properties. We 
therefore study the SSE in the ID Hubbard model, and investigate whether its spectral 
weight can be comparable to the CDE as observed in the experiments. 

In this paper, we want to study the ID single band Hubbard model (HM) through the 
Bethe-ansatz equations and the Lanczos-Gagliano method, which is shown to be in excel- 
lent agreement with the exact diagonalization result. Although HM is a lattice model of 
short-ranged on-site interaction, unlike the realistic ID QWR systems which are continuum 
systems with long-ranged Coulomb interaction, we could still use it as a valid model in 
qualitatively discussing the problem of the polarized Raman scattering because all ID inter- 
acting systems belong to the LL universality class and generic issues may be addressed in 
any particular ID model. We first identify the holon excitation of the Hubbard model to be 
the CDE, and the triplet spinon to be the SDE in the usual RRS language, by comparing 
their dispersion relations in the whole spectrum. We then further obtain the finite spectral 
weights of SSE in the charge density spectrum and show that the weight of singlet spin 
density excitation is still rather low in the HM and cannot produce large spectral weights in 
the polarized RRS scattering spectrum as found in the experiment. Thus, our HM results 
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shows that the " SPE" peak in the polarized RRS experiments is unlikely to be explained by 
the ID singlet spin excitation, at least within any nonresonant theory which considers the 
elementary excitations only in the conduction band. 



A. Theory 

The simple ID single band Hubbard model, 

H =-tH (4fv<V + H.c.) + £/5>; T n i4 , (20) 

i,cr i 

where q a and Tli (j £11*6 respectively the fermion creation operator and the density operator for 
site % and spin a, has been extensively studied, t and U are hopping energy and on-site short 
range spin- dependent interaction following the usual notation in the literature pl|-p^1 . Note 
that the Hubbard model is basically a model with a spin-dependent short range (on-site) 
interaction, U. It has generic LL properties in the long wavelength limit and for low-lying 
excitation energy, i.e. no single particle behavior in the spectral function at Fermi wavevec- 
tor. The explicitly spin-dependent interaction, U, in the HM, however, should make the spin 
singlet state more enhanced in the spectrum, and easier to study. Among the many accurate 



and useful methods to study the ID HM, we use the Bethe-ansatz method ||21| , |23| to obtain 
the ground state and the low-lying excitation state dispersion spectra. It is well-known 
that the Bethe-ansatz wavefunctions are not particularly useful in calculating correlation 
functions, and therefore we need an alternative method to obtain the spectral weights of 
the elementary excitations. We calculate the charge density and spin density correlation 
functions, to be compared respectively with the polarized and depolarized inelastic light 
scattering spectra, by using the Lanczos-Gagliano method [f2~4| -|27|j . This method (described 
and discussed below) gives a simple but fast- convergent result for the correlation functions 
in the lattice model. By comparing the momentum-energy dispersion relations of these two 
different methods, (i.e. the Bethe-ansatz and the Lanczos-Gagliano method) we can iden- 
tify each important spectral peak obtained by the Lanczos-Gagliano method to be a certain 
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Bethe-ansatz elementary excitation in the Hubbard model language (holon, triplet spin or 
singlet spin excitations, for example) and then estimate their relative weights for compari- 
son. Our results obtained by this technique are consistent with the quantum Monte Carlo 



calculations [28] where appropriate. 



1. Bethe-ansatz equations 

It is well-known that the ID Hubbard model can be solved exactly |2T| , |22[| by the Bethe- 
ansatz method. The eigenvalue equation of Eq. (20) is proved to be identical to solving the 
coupled system of equations (under periodic boundary condition) 

ik s L = tt sinfc 3 -- A Q + rt//4 

" X a -smk j +iU/4 il A a - A/3 + iU/2 

j^i \a -smkj -iU/i \ a - \p-iU/2' 1 ' 

where L(N) is the total number of sites(electrons) and M is the number of down-spin elec- 
trons (M < N/2). The pseudo-momentum, {kj}, and spin rapidities, {A a }, are generally 
complex variables to be solved and related to the physical states of energy, E, and momen- 
tum, p, by 

JV 

E = —2t ^2 cos kj, (23) 

3=1 

and 

P = J2 k 3- (24) 

3=1 

If the k/s and A a 's are all real, the identity of the phases in Eqs. (21)-(23) can be obtained 
by taking the logarithm. Then we have the following well known results, 

Lkj = 27r/, + 2 £ tan" 1 ( Xa ~jf* kj ) (25) 



a=l 

N 



U/A 



^ i / X a — sin kj \ 
2 Eton" 1 a ^ =2nJ a 



3=1 



M I 

+2 £ tan" 1 

8=1 \ 



X a — ^f: 

U/2 



(26) 
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where the quantum numbers, {Ij}, are all distinct from each other and are integers if M is 
even and are half-odd integers if M is odd, and are only defined in < L. Similarly, the 
set {J a} are all distinct and are integers if N — M is odd and half-odd integers if N — M 
is even. The value of {J a } is restricted by \ J a \ < (N — M + l)/2. Generally, it is not hard 
to use the Bethe-ansatz equations to solve large size systems. In the thermodynamic limit, 
L — > oo, one can find the equivalent linear integral equations for the density of fc's and A's 



on the real axis [P2| , P9|| . But we will only focus here on the finite size systems in order to 
compare the Bethe-ansatz results with the results of the Lanczos-Gagliano method, which 
is necessarily computationally restricted to small system sizes. 

To solve these Bethe-ansatz equations, we first have to define the proper quantum num- 
bers, {Ij} and {J a }; then solve Eqs. (25)-(26) to get kj's and A Q 's, and then get the 
momentum, p, and the energy, E, of that state specified by those quantum numbers. Here 
we present the quantum number structures of the ground state and two low lying excited 
state, the "4&f" singlet states, the "2kp" triplet states, and the "2kp" singlet states as first 
named by Schultz The first two have fc's and A's all real in Eqs. (25)- (26), while the 



last one has one pair of complex A's in Eqs. (23)-(24) 

Ground state: It is easy to see that the ground state is nondegenerate only if N is of 
the form 4m + 2 (m is an integer). In the following, we just study the nondegenerate case for 
simplicity. Considering the essential symmetries, one can write the ground state quantum 
number satisfying the above restrictions to be 

{Ij} = {-(N - l)/2, • • , (N - l)/2}, (27) 
{J a } = {-(N/2 - l)/2, • • (N/2 - l)/2}. (28) 

"4fcp" singlet state (holon excitation) : The first simplest excited states are obtained 
by removing one of the momentum quantum numbers, — (N — l)/2 + Zq, in {Ij} and adding 
a "new" one at (N — l)/2 + Iq outside the ground state sequence. All other momentum 
quantum numbers and spin quantum numbers are kept the same as in the ground state 
structure. Therefore, the new sequence of {Ij} is 
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{/.} = { _ (iV - i)/ 2 , ...,-(N- l)/2 + ^0-1, 
-(iV - l)/2 + 1 + l, . . (JV - l)/2, 

(iV-l)/2 + / }, (29) 

and the {J a } is the same as the ground state in Eq. (28). Therefore, there are two free 
parameters, io and Io, for this type of excitations. In this paper, we use (io, Io) to denote 
this excitation state. According to Schultz |14| , they are named "4k F " singlet states due 



to their energy minimum at k = 4kp in their dispersion spectrum. In the literature, these 



states are also called "particle-hole excitation" or "holon" excitation [ft2jl . In the rest of this 
paper, we will call them "holon" excitations for simplicity. (This is related to the CDE of 
our earlier sections.) 

"2kp" triplet state (triplet spinon excitation): Next we consider the excitations of 
the J's with all A's and fc's real. This is possible only if M < N/2. The simplest excitations 
of this type are obtained by considering M = N/2 — 1. The total spin of the system is 
S — 1, so we expect this excitation to be related to a triplet spin excitation. The quantum 
numbers of these states are 

{I j } = {-N/2 + l,---,-N/2 + i -l, 

-N/2 + i + l,--;N/2,N/2 + I }, 
Ji = -N/4 + 6 au i, 

where a = 2, M and 1 < «i < «2 ^ M + 2. Here ot\ and «2 are the free parameters in the 
spin quantum number, { J a }, and io and Jo are the two parameters in momentum quantum 
number, {Ij}- From Eqs. (23)-(26), we can see that io and Iq shift the total momentum and 
energy of the spectrum created by the spin excitation in { J a }. In the following calculation, 
we use (io, Io, a\, 0:2) to denote these states in the spectrum. These excitations are called 
"2kp" triplet states because its minimum energy is at k = 2kp. In the rest of this paper, 
we will call them "triplet spinon" for simplicity. (This is related to the SDE of the earlier 
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sections.) 

"2A; F " singlet state (singlet spinon excitation). The third possible elementary ex- 
citations are from the complex solutions of Bethe-ansatz equations, Eqs. (21)-(22). Spin 
singlet states (M = N/2 and then S = 0) are obtained by having one pair of the complex 
conjugate, A± = A_r ± A/ with all other /c's and A's real. The new set of Bethe-ansatz 
equations are obtained to be 



where j = 1, 2, N and a = 1, 2, M, but a ^ a±, a 2 - The two equations for the complex 
A± are 




(31) 




(32) 




( 



(XR-smktf + jU/A + Xj) 2 
(XR-smk.y + iU/A-Xj) 2 



) 




(33) 



= 2ttJ 




(34) 



where 
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integer < 



if |A/| > £7/4, and iV - M is even 

or if | A/ 1 < 17/4, and N — M is odd (35) 
half odd integer otherwise. 

As for the quantum number, {Ij} and {J a }, in the singlet states, we choose them to be 
the same as the ground state, Eqs. (27)- (28), except for the two free "holes" at J Q1 and 
J a2 , whose related spin quantum numbers, A Ql and \ a2 are replaced by the pair of complex 
conjugate, A± = A#±Aj. Eqs. (33)-(34) are usually too complex to give a nontrivial solution 
because the usual numerical iteration method will converge to the trivial A/ = solution. 
But one can simplify these equations by taking A/ = U/4 and one kj satisfying sin kj = \r, so 
that Eq. (33) could be neglected, and all the terms containing tan _1 [(A# — sin kj)/(U/A — A/)] 
in Eq. (34) contribute a phase ±tc. The phase number, J, is set to make the total phase 
shift (including those from tan _1 [(A_R — sin kj)/{U/A — A/)]) to be zero in the calculation. 

The spin singlet excitations have a dispersion similar to the triplet ones. Here we could 
use (0:1,0:2) as the quantum number to define these states. In the finite size system with 
repulsive interaction, U, the singlet states have higher energy than the triplet ones, but 
they will become degenerate in energy as we go to the thermodynamic limit (L — > 00, and 
(n)=constant). In the experiment, the spin triplet excitations (i.e. SDE) are observed in the 
depolarized RRS spectra where a net spin-flip occurs, while the singlet states are observed 
in the polarized spectra, which involve no net spin-flip. 

As mentioned in the beginning of this section, the Bethe-ansatz method does not, in 
general, provide the spectral weights for their solutions. Therefore, the three elementary 
excitations above may not be equally important from the experimental point of view, i.e. 
they may carry very different spectral weights (and some may even be unobservable in the 
experimental spectra). All we know from the Bethe-ansatz solutions are the existence and 
the dispersion of these excitations but not their spectral weights. We also know that these are 
allowed solutions of the HM just as the SPE is an allowed solution of the FL-RPA model (but 
not the LL model). Comparing the mode spectral weights calculated by Lanczos-Gagliano 
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method we discuss next, we calculate the relative spectral weights of these solutions and 
then study their interaction dependence. 



2. Lanczos-Gagliano method 

In this paper Lanczos-Gagliano method means the combination of two important tech- 
niques in the lattice model. The standard Lanczos method is to construct an L x L matrix 
representation for the tridiagonal Hamiltonian, like Eq. (20), and then directly diagonalize 
it to get the eigenvalues, E n , and eigenfunctions, <& n , which could be used to do further 
calculations, such as obtaining spectral weights. But since only ground state energy and 
wave function are needed in calculating the correlation function by using the Gagliano's 
method (see below), we use a simpler but more efficient way, the modified Lanczos method, 
to calculate the ground state energy and wave function. This method has been analyzed 
and discussed in detail in references pl-|2"7J|, and here we only provide a brief review for the 



sake of completeness. 

In the modified Lanczos method, we first present H in a 2 x 2 matrix in the basis of ipQ 
and ipii which are both normalized and orthogonal with each other by setting 

h = ((g 2 > _ (g>2) i/ 2 - ( 36 ) 

where ipo could be an arbitrary trial wave function, and (H n ) = (tpQ\H n \ifjQ) . After diago- 
nalizing the 2x2 matrix representation of H, we can get an eigenenergy, £q, and a wave 
function, ipo (subscript means the lower energy eigenstate). This Eo and "00 are then a 
better approximation to E and $o than the initial trial state, (H) and ipQ. Using the same 
process by iterationally replacing i/jq by ipo to construct a new 2x2 matrix representation 
for H, we will then get a more accurate approximation to the ground state energy and wave 
function simply by diagonalizing the 2x2 matrix. 

The Gagliano's method p4Vf2"T| is to use an infinite continued fraction representation to 



calculate a general dynamic correlation function, XaOZj^Oj where A could be any operator. 
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Its advantage is that one just needs the ground state information to get the dynamical 
results extremely efficiently in a very simple and straightforward way. The convergence and 
accuracy of the calculation is very good compared with the much more time consuming exact 



diagonalization method [Z6[|. In this paper, we want to study the charge density correlation 
function, XpiQ,^), an d the spin density correlation function, x CT (g, cu). Their corresponding 
operators, A's, are respectively defined as 

Mm) = P(Q) = 4ff E e- lQl feW - (n>) , (37) 



and 

A(q) = a(q) = -±= £ e^Vc^c^, (38) 

where (n) = N/L is the average density. The Gagliano's method introduces a function, 
G A (Z), in the following formula to calculate the correlation function 

XA (q^) = ^(* \A\q)\* n )(* n \A{q)\* ) 

n 

x5(u-(E n -E Q )) (39) 
= -ImG A (u + Eo + irf, (40) 

7T 

where 7 is a (small) phenomenological broadening factor, and the function, Ga is given by 



G A (Z) = {%\A\q){Z-H)- l A{ q )\%) 
(^ \A^q)A( q )\%) 

z ~ a ° W~ 



(41) 



Z — ai 



2 



Z 

The coefficients, {a n } and {b n } are defined by generating a set of orthogonal states with the 
following iteration formulae, 

\f ) = A(q)\<t> ) 

\f n+1 ) = H\f n )-a n \f n )-b 2 n \f n ^) (42) 

a n = (fn\H\f n )/(f n \f n ) (43) 
b 2 n+1 = (fn + i\f n+ i)/(fn\fn), and b = . (44) 
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Therefore, after we get Eq and $o from the modified Lanczos method, we can calculate {a n } 
and {b n } iteratively and then obtain the correlation function xpiQ^) an d X<t(<?, ^) by using 
Eqs. (39)- (44). Our calculation of the expectation value is done in the zero temperature 
ground state of the system. The method may be easily generalized to finite temperatures [|27| . 
which is, however, unnecessary for our purpose of studying the ID excitation spectra. 



B. Results and discussion 

We study the ID Hubbard chain with three different densities, (n) = N/L = 1/3 for 6 
electrons in 18 sites, (n) = 1/2 for 6 electrons in 12 sites, and (n) = 5/6 for 10 electrons 
in 12 sites. (Note that the usual filling factor of the system is (n) jl since our definition 
of density does not include spin.) The size of the Hubbard chain is dictated here entirely 
by computer memory restrictions in calculating the correlation function via the Lanczos- 
Gagliano method. We keep the electron number to be 4m + 2 with m an integer in order 
to have a nondegenerate and zero-momentum ground state under the periodic boundary 
condition. Throughout our calculations, we set the broadening factor, 7 in Eq. (40), to be 
O.Olt (where t is the nearest-neighbor hopping amplitude in Eq. (20)) and use the modified 
Lanczos method to calculate the ground state energy iterationally until convergence to within 
less than 0.1% in the ground state energy is reached. We also truncate the infinite continuous 
fraction in Eq. (41) at 25-27th order terms, which gives us good convergent results in the 
calculation. 

In the following, we will first discuss the results related to the polarized spectrum, which 
involves no net spin-flip in the system, by using the two methods mentioned above at a 
fixed interaction strength, U/t = 3. Then we consider the depolarized spectrum under the 
same conditions. Finally we discuss their interaction, U, dependence by varying U/t in our 
calculations. In the discussion below the terms "resonance" or "resonance peaks" refer to 
the Lanczos- Gagliano calculations. 
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1. Polarized spectrum analysis 



We compare the dispersion of the charge density excitation with the dispersions of the 
"4/cp" singlet states (holon) and the "2kp" singlet states (singlet spinon) given by the solu- 
tions of the Bethe-ansatz equations, because these two are the low lying elementary spinless 
(S = 0) excitations of the ID Hubbard model and as such should correspond to the polarized 
spectrum. We will also study their relative spectral weights. Both lower density ((n) = 1/3) 
and higher density ((n) = 5/6) results are shown together (Fig. 6-9) for further discussion. 

In Fig. 6(a), we show the spectral dispersion obtained by the poles of the imaginary 
part of the charge density correlation function. The doping density is (n) = 1/3 for 6 
electrons in 18 sites in the ID Hubbard chain with periodic boundary conditions. The 
center of each open diamond represents the position of the pole, and its area is proportional 
to the spectral weight of that excitation. In the same figure, the dispersions of the holon 
(star) and singlet spinon (open square) excitations given by the solutions of Bethe-ansatz 
equations are also shown for comparison. Several features could be noted from Fig. 6: 
(i) the excitation peaks of the charge density correlation function have a linear dispersion 
in the long wavelength limit (q <C kp = 7r/6) and its slope gives the velocity of charge 
density excitation of ID Hubbard model, (ii) The resonance peaks form a wing upto the 
large momentum region (i.e. low energy excitations correspond to the low momentum and 
high ones to high momentum), with a maximum energy uo = At at q = tt. (iii) The sizes 
of the diamonds, which represent their spectral weights, show that the peaks at higher 
energy generally have greater spectral weights than the ones at lower energies at the same 
momentum (i.e. spectral weights are greater for peaks at higher energies). Therefore one 
could see by eye a sine-like curve at the upper edge of the resonance wing with a maximum 
at uo = At, and this observation is consistent with the results from quantum Monte Carlo 



simulations on larger system [p8~ |. (iv) There are no excitation states for singlet spinons at 
small g = 7r/9. This implies that the singlet spinon of the ID Hubbard model is not allowed 
for momentum smaller than 2 x 2tt/L, where 27r/L is the momentum scale of this finite size 
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(L) system. This follows from the fact that the singlet state must be excited by a pair of 
complex conjugate A± in Eqs. (21)-(22), which is at least a two-particle excitation, so that 
the minimum momentum required is 2 x 2rr/L. (v) One interesting feature is that there 
are clearly two energy minima at q — 2/cp = tt/3 and q = 4kp = 27r/3 in the spectrum. 
Comparing these resonance peaks with the solutions given by Bethe-ansatz equations, we 
find that the holon excitations cover almost exactly the same region including the energy 
minimum at 4kp except for the lower-lying peaks around 2kp, where the singlet spinon just 
matches those peaks. In other words, we could reasonably claim that the most dominant 
features of the resonance peaks given by the charge density correlation function arise from a 
combination of holon and singlet spinon excitations in the ID Hubbard model. This result 
could not be trivially obtained either by solving Bethe-ansatz equations or by calculating 
the charge density correlation function alone, as mentioned in the introduction — one must 
combine the two techniques to come to this conclusion. Other spin singlet excitations given 
by the solutions of Bethe-ansatz equations (for example, 2 pairs of complex A's in Eqs. (21)- 
(22)), carry very small spectral weights because no other significant resonance peaks are 
found in this dispersion spectra, except for some trivial ones. In the thermodynamic limit, 
we expect that only the 4kp holon and 2k p singlet spinon will have finite spectral weights 
and could be interpreted as the "charge density excitation" and "single particle excitation" 
in the RRS spectra respectively when comparing with the experiments ||| as we mentioned 
in the earlier sections. We discuss this issue further later in this paper. 

In Fig. 6(b), we show the imaginary part of the charge density correlation function of 
the same system at q = 2tt/9. It shows that singlet spinons have a relatively small, but 
non- negligible weight, compared with the weight of the dominant holon excitations. Their 
relative spectral weight ratio (singlet spinon/holon) is less than 0.1. Similar results are also 
obtained in the systems of 6 electrons in 12 sites, (n) = 1/2, which are not shown in this 
paper. 

In Fig. 7(a), we show the dispersion of the resonant poles of the charge density correlation 
function of 10 electrons in 12 sites ((n) = 5/6). Here the holon excitations form a more 
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narrow wing than in the lower density system, but the basic shape of the dominant curve is 
almost the same. Below this curve, the singlet spinon occupies almost the whole resonance 
region. Since 4fcp = 5n/3 in this high density system, we cannot see the 4k f energy minimum 
in this figure (actually, the gap of this energy minimum is very large in this finite system, 
but will go to zero in the thermodynamic limit [|14|| ). But one could still see the energy 
minimum of the singlet spinon at 2kp = 5ir/6 ~ 0.8337T in Fig. 7(a). There is a notable 
feature in these results at energies above the holon excitations: there are some additional 
states, having an energy gap equal to U(= 3t) and a maximum energy greater than 4t at 
q = 7r in this high density situation. These states (solitons) must arise from the double 
occupancy of the electrons in the Hubbard model, which consequently explains their high 
energy status. Quantum Monte Carlo method, to the best of our knowledge, does not 
provide any information about these high energy double occupancy states at the same density 



n = 5/6 in the literature ||28|| . From the Bethe-ansatz equations point of view, however, 
these states should be obtained by taking the complex momentum, kj, solutions in Eqs. 
(21)-(22) whether in the high or the low density system. Once again, we see the 
importance of studying the spectral weights of the Bethe-ansatz solutions by comparing them 
to the correlation function results so that one could tell the most realistic and physically 
meaningful states. Just having the solutions, without much idea about their spectral weights, 
is not useful in determining the experimental and/or physical relevance of the particular 
excitations. Therefore, as shown in Fig. 7(b), the three most important contributions to the 
charge density correlator arise from singlet spinon, holon, and soliton excitations (double 
occupancy excitations), from lower energy to higher energy regime respectively in the ID 
Hubbard chain. Their relative spectral weights show that the singlet spinon has the smallest 
weight, and it could be shown that there are no gapless holon and singlet spinon excitations 
in the half-filling ((n) = 1) ID Hubbard model systems, where the soliton (double occupancy 
excitations) and other higher energy states dominate the excitation spectrum. 
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2. Depolarized spectrum analysis 



In Fig. 8(a), we show the resonance dispersion spectrum of the spin density correlation 
function, (era), of the low density system (6 electrons in 18 sites). The triplet spinon 
excitation spectrum given by the solutions of Bethe-ansatz equations is also presented for 
comparison. Several features are found: (i) in the long wavelength limit, the resonant poles 
have a linear dispersion, whose slope gives the velocity of the spin density excitation. One 
could easily see that this velocity is always smaller than the velocity of the charge density 
excitations at the same density. This is consistent with result of previous work (ii) The 



resonance poles form a wing upto the large momentum region (q ~ it), whose maximum 
excitation energy is below At. (iii) Unlike the results for the charge density correlation 
function in Fig. 8(a), the most dominant poles are located in the lower energy part of the 
resonance wing, which correspond to the triplet state without any excitations in {Ij}, and 
therefore is related to the lowest energy ones in our calculation, (iv) The resonance spectrum 
has an energy minimum at 2kp. Compared to the triplet solutions of Bethe-ansatz equations, 
the triplet spinon excitation spectrum has only three peaks matching the resonance poles of 
the largest spectral weight at their momentum values (q = tt/9,2tt/9, and 3tt/9), and the 
other three match the poles of relatively much weaker states (not visible in Fig. 8(a), but 
distinguishable in the absorption spectrum shown in Fig. 8(b)). This result demonstrates 
that the spectral weights of elementary excitations could be very different even if they result 
from the same type of the Bethe-ansatz solution. Fig. 9 shows the dispersion relation of 
the spin density correlator of the large density system ((n) = 5/6) and the corresponding 
triplet spinon excitations by the Bethe-ansatz solutions are also shown. 

3. Interaction dependence 

In this section, focusing on the lower density system (L = 18 and iV = 6) and a fixed 
momentum (k = 2tt/9), we study the mode dispersion and spectral weight of these excita- 
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tions in a range of finite interaction (U/t < 10) to obtain the interaction dependence of the 
excitation spectra. First, we study the polarized spectrum given by the imaginary part of 
the charge density correlation function, (pp) (shown in Fig. 10(a)). Then we compare the 
energy of the resonance peaks in the series of spectra with the Bethe-ansatz results (Fig. 
10(b)), and discuss the interaction dependence of the mode velocity (Fig. 10(c)). Finally 
we discuss the interaction dependence of the spectral weight for each elementary excitation 
(Fig. 11). 

In Fig. 10(a), there are basically three peaks in the typical structure of the polarized 
spectrum, and we can identify them as the singlet spinon, the second and the first holon 
excitation (from lower to higher energy) by explicitly comparing with the energy given by the 
Bethe-ansatz solution in Fig. 6(a). Using the notation introduced in Sec. IV- Al, the singlet 
spinon is the state (2,3), while the holon I and II states are (6,2) and (5,1) respectively. 
Several interesting features can be found in Fig. 10(a): (i) in the noninteracting (U = 0) 
case, there are only two equal weight poles, which could be understood as the two single 
particle (electron and hole) excitations around Fermi surface, k = k F . (ii) When finite 
interaction is turned on, there is an additional excitation. According to the comparison of 
dispersion relations, both the new peak and the higher energy peak should be interpreted as 
holon excitations (called holon II and holon I respectively, corresponding to different {ij}'s). 
(iii) The singlet spinon excitation (shown in Fig. 10(a)) has a rapidly decreasing spectral 
weight with increasing interaction, and disappears totally as U/t > 5.0. (iv) The two holon 
excitations shift to higher energy as U increases, and maintain almost the same spectral 
weight except one more peak appears as U/t > 8.0 (see Fig. 10(a)). Above U/t = 8.0, the 
appearance of the new small peak affects both the spectral weight and excitation energy 
of the holon II excitations (see Figs. 10(b) and 11(a)). There are basically two possible 
interpretations for this result. One is that this peak does not represent real excitations, 
but may be arising spuriously from the inaccuracy of the finite truncated series or finite 
iteration used in the Lanczos-Gagliano method in the large U range (see Eq. (41) and the 
discussion in the last section). Another possible reason is that it may arise from the higher 
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energy excitation states of unknown origin, which are also obtainable from the Bethe-ansatz 
solutions, but whose strength is visible only when the interaction strength is large enough. 
We will not discuss this anomalous peak any further in this paper since this falls outside the 
scope of our main interest. 

In Fig. 10(c), we plot the excitation velocity, which is defined to be 

as a function of interaction strength. We find that when the interaction is weak (U/t < 1), 
the two (holon and singlet spinon) excitations are almost degenerate, while their relative 
spectral weights change a lot (see Fig. 11(a)) as a function of U/t. When U/t increases, the 
holon has greater excitation energy and hence velocity, but the velocity of the singlet spinon 
decreases fast. This result holds even in the thermodynamic limit. 

In Fig. 10(b), we see more clearly that the energies of the three elementary excitations 
are only weakly dependent on the interaction for small U, but strongly dependent on U 
for large U. In Fig. 11(b), we have a logarithmic scale in the spectral weight dependence 
with respect to the interaction in small U/t range (U/t < 1.0). By calculating the slope of 
these data, we find that the spectral weight of holon I, Shoion i, is almost a constant in the 
smaller interaction range (U < 0.5t), and then weakly decreases for higher U. However, the 
spectral weight of holon II has a stronger power-law behavior, Shoion n oc U imb . Thus the 
two holon excitations differ a great deal in their interaction dependence of their respective 
spectral weights. 

In Fig. 12(a) we show the calculated depolarized spectra by taking the imaginary part 
of spin density correlation function, (cxcx) for various interaction strengths. In the nonin- 
teracting case, the spectrum is the same as the polarized one in Fig. 6(a) due to spin 
rotational symmetry. But with increasing interaction strength both triplet spinon peaks 
move to lower energy in contrast to the polarized spectra. Compared with the Bethe-ansatz 
results in Fig. 12(b), the lower/higher energy peak, triplet spinon I/II, is the state de- 
noted by (6,0,1,3)/(6,1,2,3). The excitation energy of the triplet spinon we obtain by the 
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Lanczos-Gagliano method agrees well with the Bethe-ansatz result in general except that 
the energy of the triplet spinon II does not seem to agree well when the interaction is larger 
than U/t = 3. From the result in Fig. 12(a), we can see that this may be due to the appear- 
ance of another excitation peak in the Lanczos-Gagliano spectra, which is not represented 
in our Bethe-ansatz solutions. Based on these results we conclude that the triplet spinon 
excitations are likely to be the dominant contributions in the depolarized spectra. In Fig. 
13, we show the interaction dependence of the spectral weights of the two triplet spinon ex- 
citations. The triplet spinon I(II) has a maximum(minimum) spectral weight at some finite 
interaction, U/t — 3 ~ 4, and the interaction dependence of the spectral weight is nontrivial. 
This result demonstrates the importance of the intermediate interaction strength of the ID 
Hubbard model. 



C. Discussion 

We systematically study the elementary excitations of ID Hubbard model by combin- 
ing the techniques of the exact Bethe-ansatz equations for the mode dispersion and the 
Lanczos-Gagliano method based spectral weight calculation of the correlation functions. 
Three types of elementary excitations, holon, singlet spinon, and triplet spinon excitations, 
are studied at zero temperature and different densities ((n)=5/6, 1/2 and 1/3) and differ- 
ent interaction strength, U. We first compare the energy-momentum dispersion relations of 
these excitations obtained by both methods and then study the mode spectral weights in 
different situations. The comparison between Bethe-ansatz solutions and resonance peaks 
of the Lanczos-Gagliano correlation function gives us some important results: (i) the holon 
and the singlet spinon excitation states show up together in the charge density correlation 
spectra. Holon states have higher energy with an energy minimum at k — Akp while the 
singlet spinons lie in the lower energy region with an energy minimum at k — Ikp. There are 
no other states of prominent spectral weights except the gapped double occupancy (soliton) 
states near half filling. This result connects the theoretical calculation of the ID Hubbard 
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model with observable physical quantities — in particular, these are the only two modes 
which are likely to show up in the polarized Raman scattering experiment probing the charge 
density excitation spectra. Another implication of this result is that one can interpret the 
Bethe-ansatz quantum numbers, {Ij} and { J a }, as the ones of collective excitations. But our 
spectral weight analysis shows that most of the Bethe ansatz solutions for the ID Hubbard 
chain do not have any observable contributions to the real physical quantities because they 
carry essentially no spectral weights, (ii) The excitation holon II has a power-law behavior 
in its spectral weight with respect to the interaction strength in the small U/t region, while 
the holon I has almost interaction- independent spectral weight (here the holon I/II could 
be generalized to represent the 4/c^-singlet excitations having holes in the edge/middle of 
the charge quantum number, {Ij})- An interesting problem for further research is to ob- 
tain an analytic formula for the exponent of the holon II excitation. This will relate to 
the small interaction expansion of Bethe ansatz equations and wave functions, which have 
not yet been explored much in the literature. When the interaction strength increases, on 
the other hand, the spectral weights of these two holons become equal as shown in Fig. 
10(a). (iii) As for the singlet spinons, we find that their spectral weights decrease to zero 
very fast (exponentially) as U increases. This could be understood from the fact that the 
on-site repulsive interaction, U, prevents the formation of the symmetric electron orbital 
wave functions, which must accompany the antisymmetric spin singlet wave function, thus 
suppressing the singlet spectral weight for large U . (iv) From the imaginary part of the 
spin density correlation functions we find that the triplet spinon is the only low energy spin 
excitation in the long wavelength limit. There is no other excitation of important spectral 
weights in this region, (v) The spectral weight study shows that among many triplet spinon 
excitation states, only those with some some special quantum numbers could possibly have 
relatively greater weights at a given momentum (see Fig. 8(a) and the text related to that) 
for finite interaction strength, U/t. Others have very small or trivial weights, which are not 
physically significant, (vi) Finally the interaction dependence of the spectral weights of the 
triplet spinon I and II differs very much in the intermediate interaction range, but becomes 
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similar in magnitude in both the weakly interacting and the strongly interacting situations. 
This shows the subtle complications in interpreting various excitation modes in the ID Hub- 
bard model for intermediate interaction strength (say U/t ~ 3). Further research is needed 
to provide a more complete understanding of this intermediate interaction region, and our 
results should be considered a preliminary investigation. 

V. CONCLUSION 

In this paper, we calculate the ID charge and spin density correlation functions (Eqs. 
(3) and (4)), which respectively correspond to the polarized and the depolarized Raman 
scattering spectra, in three different models: Fermi liquid-RPA model, the Luttinger liquid 
model, and the Hubbard model. Fermi liquid results (in RPA) show a strong CDE mode 
(plasmon) with a very weak SPE mode in the polarized spectra, while the Luttinger liquid 
results show the bosonic collective excitation (CDE or plasmon) as the only mode in the 
polarized spectra without any SPE. The ID Hubbard model results show that both holon 
and singlet spinon excitations contribute to the polarized spectra, and we study their relative 
spectral weights and interaction dependence in detail. Note that the Hubbard model being a 
lattice model, its applicability to the RRS experiments in ID semiconductor quantum wires 
(which are continuum systems) is necessarily limited and is at best qualitative. 

Focusing on the polarized RRS calculations, we can compare the lattice HM results with 
the continuum FL-RPA and LL results in the context of experimental observations in ID 
Q WR systems [[§] . According to the solutions of the HM Bethe-ansatz equations, we find 
that the RRS spectrum is best interpreted as the combination of holon and singlet spinon 
excitations as well as the CDE and "SPE" (or "SSE", i.e. the spin singlet excitations) in the 
high and low energy regions respectively. Therefore, the numerical ratio of spectral weights of 
the singlet spinons to the holons should be interpreted as the ratio of the SSE(SPE) spectral 
weight with respect to the CDE spectral weight. The experimental results show that SPE 
has almost the same spectral weight as the CDE, i.e. this ratio is almost unity, whereas 
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in our calculation, this ratio is much smaller than unity in the intermediate interaction 
range (U/t > 0.5). For weak interactions (U/t < 0.5) it seems from Fig. 11 that SSE 
could have comparable weight to CDE. But this raises a serious problem with respect to the 
experimental results Q , which find the plasmon velocity to be at least two times greater than 
the Fermi velocity in the experimental long wavelength regime, while this is not true for the 
ID Hubbard chain unless U/t > 1 as seen from Fig. 10(c). Thus the plasmon dispersion (if 
we demand agreement between HM and the experimental results) indicates the appropriate 
interaction range in the HM to be U/t > 1, where our calculated SSE/CDE spectral weight 
ratio is extremely small, in contrast to the experimental finding. Therefore, we can conclude 
that the ID Hubbard model results can not explain the experimental data consistently for 
both the mode dispersion and the spectral weight, although the lattice nature of the Hubbard 
model makes a comparison with QWR experiments somewhat misleading. Moreover we note 
that in the Hubbard model, the on-site interaction, U, is short-ranged (in fact, zero range), 
and operates only between opposite spins at the same site. We therefore expect that the 
spectral weight of the singlet spinon is (perhaps strongly) overestimated in the Hubbard 
model compared with the Coulomb interaction (which is explicitly spin- independent) case. 
We believe the spectral weight of SSE (or SPE) in a continuum interacting model should 
be even less than that in the finite-sized Hubbard model here, just as we find in the RPA 
(or the HA) results in the FL theory where the mode dispersion (both for CDE and SPE) 
agrees well with experiments, but the spectral weight for the SPE is far too weak in the 
theory. Based on this logic our results seem to rule out the interpretation that the SSE 
shows up as the additional anomalous SPE mode in the RRS spectrum. The singlet spin 
excitation according to our analysis could not provide sufficiently large spectral weights in 
order to quantitatively explain the large spectral weight in the so-called SPE mode seen in 
the experimental RRS spectroscopy of ID QWR systems. This conclusion, as emphasized 
before, only applies to the NRS nature of our theory — if the valence band and the resonant 
nature of RRS experiments is playing an important role, as likely J^-fl? then the theory 
needs to take resonance into account. 
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In summary, we systematically calculate the charge density (polarized spectra) and the 
spin density correlation (depolarized spectra) functions of one-dimensional systems in three 
different models: Fermi liquid model, Luttinger liquid model, and Hubbard model. In 
the polarized spectra, we find that the FL model shows a strong collective charge density 
excitations at plasmon energy and a relatively weak single particle excitation at u = qvp, 
while the LL model shows one bosonic (plasmon/CDE) excitation only. Comparing the 
plasmon excitation energy of FL model and the bosonic excitation of LL model (see Eqs. 
(11) and (17)) we find these two excitations are identical, and the small SPE peak in FL 
model arises from the finite curvature effect of electron energy dispersion at the Fermi 
point. In the Hubbard model, however, two excitations, holons and singlet spinons, show up 
together in the polarized spectra. We show that the holon excitations are actually the CDE 
in FL model or the bosonic excitation in LL model, while the singlet spinons in the HM 
arise from the spin degree of freedom and finite dispersion curvature at Fermi point. If we 
compare the spectral weights of the lower energy excitations (SPE of FL model/no peak in 
LL model/singlet spinons in Hubbard model) and the weights of the higher energy excitations 
(CDE in FL model/boson peak in LL model/holons in Hubbard model), we find that the 
higher energy excitations always have (much) larger spectral weight than the lower energy 
ones in all models. This shows that the equal weight two-peak structure observed in the 
experiments Q could not be explained by the nonresonant Raman scattering mechanism, no 
matter how one interprets the lower energy excitations to be SPE or SSE. Recent theoretical 
work |7]-|9| on resonant Raman scattering spectroscopy indicates that the low energy SPE 
feature may be a purely band structure effect arising from the participation of the valence 
band in the resonant scattering process. This also explains why this anomalous peak shows 
up in all dimensions in experiments and not just in ID. In the depolarized spectra, however, 
only one spin excitation (the SDE or the spin triplet excitation) is observed in these three 
models. The vertex correction of the FL model will in general reduce the SDE energy 
compared with the SPE energy, and separate the SDE from the SPE. In the intermediate 
interaction region, the two triplet spinons in the Hubbard model have very different spectral 
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weight behavior, showing very interesting interaction effects which need to be studied in 
more details in the future. 
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FIGURES 




FIG. 1. Diagrammatic representation of the conduction band irreducible polarizability, 
Ho(q,u) and reducible polarizability, H(q,u>), in standard RPA calculation. V c (q) is the Coulomb 
interaction. 
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FIG. 2. (a) The energy-momentum dispersion relation for the plasmon mode and the SPE 
region of ID system, (b) The dynamical structure factor of the polarized RRS spectrum in RPA 
calculation for the ID quantum wire system at q = O.lkp. Vertex correction in Hubbard approx- 
imation is also shown for comparison. Parameters are the same as the experiments in ref. [4]. 
Finite broadening factor is involved to present the delta-function peaks. 
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FIG. 3. The dynamical structure factor of the depolarized RRS spectrum calculated from 
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FIG. 4. The dynamical structure factor of the polarized RRS spectrum calculated by including 
(a) the breakdown of momentum conservation and (b) the nonparabolic energy dispersion. The 
dot, dashed, and solid lines in (a) represent the broadening parameter T = 0, 10~ 3 , and 4 x 10~ 3 
respectively (see Eq. (12)), and in (b) represent the nonparabolicity parameter A = 0, 0.02, and 
0.1 respectively (see Eq. (13)). All these effects cannot enhance the SPE spectral weight to be 
comparable to CDE in the reasonable range of V or A. 




FIG. 5. Diagrammatic representation of the Ward identity for the vertex function r rs . The 
solid lines represent the single-particle Green's function while the wave lines represent the interac- 
tion. 
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FIG. 6. (a) Energy-momentum dispersion relation and (b) the spectrum of charge density 
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portional to the spectral weight of each charge(spin) excitation peak. The numbers above the 
holon and singlet spinon peaks in (b) are the quantum numbers defined in Sec. IV-A1 from the 
Bethe-ansatz equations. 
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(n) = 6/18 = 1/3, at k = 2w/9, and (b) the excitation energies of the three elementary exci- 
tations in various U/t, compared with the Bethe-ansatz results (solid lines), (c) The velocities of 
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FIG. 11. The spectral weights of the three excitations in the polarized spectrum as a function 
of the interaction strength, U/t, in (a) linear scale from U/t = to U/t = 10, and (b) linear-log 
scale for small U (U/t < 1.0). The momentum k = 2ir/9 is the same as in Fig. 6. The inset of 
(b) is the log-log plot for the holon II excitation. We can see that the holon II excitation increases 
power-lawly in its strength as U increases, showing a possible Luttinger liquid behavior in the 
weakly interacting system (see text). 
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FIG. 13. The spectral weights of the triplet spinon excitations as a function of interaction, U/t, 
in linear scale from U/t = to U/t = 10. We consider the same system and momentum as in Fig. 
12. 
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